Collective cavity QED with multiple atomic levels 
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We study the transmission spectra of ultracold Rubidium atoms coupled to a high-finesse optical 
cavity. Under weak probing with vr-polarized light, the linear response of the system is that of a 
collective spin with multiple levels coupled to a single mode of the cavity. By varying the atom 
number, we change the collective coupling of the system. We observe the change in transmission 
spectra when going from a regime where the collective coupling is much smaller than the separation 
of the atomic levels to a regime where both are of comparable size. The observations are in good 
agreement with a reduced model we developed for our system. 
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I. INTRODUCTION 

The coherent interaction between light and matter in 
the context of cavity quantum electrodynamics (cavity 
QED) has been under intensive study in recent years, 
especially in the context of quantum information pro- 
cessing [Tj . The essence of this interaction is captured by 
considering a two-level atom coupled to a single mode of 
the electromagnetic field inside a cavity [2 [3]. However, 
for many physical systems, two important extensions to 
this model have to be taken into account. First, real scat- 
terers often have a more complex level structure [4] and 
second, the coupling for N scatterers coupled to the same 
cavity mode is collectively enhanced [5 . Both extensions 
provide useful functionality in the context of quantum 
information processing. The collective enhancement has 
recently been used to store multiple microwave pulses 
in the collective modes of an electron spin ensemble [6], 
whereas different levels could be used to collectively en- 
code several qubits [7], resulting in quantum repeaters 
able to perform simple error correction [8 . 

Studying a system of N alkali-metal atoms coupled 
to a cavity allows one to study both effects in detail. 
When probing the system with 7r-polarized light for ex- 
ample, the transmission spectrum deviates strongly from 
the simple two-level picture. For a two-level atom, the 
spectrum shows a single splitting due to an avoided cross- 
ing between the energies of the state with one excitation 
in the atom and the state with one excitation in the cavity 
[9^. Alkali-metal atoms have multiple hyperfine ground 
and excited states. In the spectrum, there are avoided 
crossings associated with the transitions between the dif- 
ferent ground and excited states [4 and the net spectrum 
depends on the relative strength of the transitions and 
their separation. For N atoms the shape of the spectrum 
depends on the atom number N as well, because the size 
of the splittings that correspond to the different avoided 
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FIG. 1. Setup of the experiment, (a) We trap an unpolarized 
gas of ultracold ^^Rb atoms in a two-dimensional lattice in- 
side the mode of a high-finesse optical cavity and probe the 
system weakly with vr-polarized light, (b) Under weak driving 
the system behaves as a collective spin with three transitions 
coupled to a single mode of the cavity. Of interest are the 
transitions to the excited states \F' — 1,2,3) of the D2 line. 
The corresponding splittings are separated by the hyperfine 
splitting of the excited states. Their size is given by the col- 
lective coupling, which in turn depends on the atom number 
N and an effective single atom coupling constant gpf. By 
varying N, we change the collective coupling and are able to 
explore different regimes of cavity QED with multiple atomic 
levels. 
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crossings is collectively enhanced by a factor vA^, as has 
been observed for thermal atoms [lQtil3j and for Bose- 
Einstein condensates [14, 15 . 

Here, we experimentally investigate the transmission 
spectra of N ^^Rb atoms coupled to a high-finesse op- 
tical cavity when the system is weakly probed with tt- 
polarized light. By changing the atom number, we are 
able to go from a regime where the size of the splittings is 
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FIG. 2. Theoretical transmission spectra for weak driving. 
Three avoided crossings between the energy of the bare cav- 
ity (soUd diagonal line) and the energies of the bare atom 
(solid horizontal lines) contribute to the transmission spec- 
trum. Transmission through the cavity is expected at the 
eigenenergies of the system. The amount of transmission is 
proportional to (a^a) of the corresponding eigenstate, as indi- 
cated by color. For low atom numbers the spectrum consists 
of three separate splittings, while for higher atom numbers 
these merge into a single splitting. Calculations were per- 
formed for an equal distributions of states and a ^ of 
6.5 MHz. All detunings are with respect to the \F = 2) to 
\F' = 3) transition. 



much smaller than the separation of the hyperfine excited 
states to a regime where both are of comparable size. 
Depending on the regime, we observe markedly different 
spectra. For low atom numbers, the spectrum shows sep- 
arate splittings, whereas for higher atom numbers these 
merge into a single splitting. In the rest of this arti- 
cle, we first extend previous models of cavity QED with 
multiple hyperfine levels [4| ITBHTS] to derive an effective 
Hamiltonian describing our system and then detail the 
experiments performed and compare them with the pre- 
dictions from theory. 



II. THEORY 

Our setup is sketched in Fig. [l] (a). We trap N ^^Rh 
atoms inside the mode of a high-finesse optical cavity. 
Our cavity is a Fabry-Perot resonator with no discernible 
birefringence. Hence, it supports two orthogonal and de- 
generate modes of the electromagnetic field. By applying 
a magnetic field perpendicular to the optical axis of the 
cavity, we ensure that these modes correspond to tt and 
±= {a-^-\-a- ) / transitions of the atom. The magnetic 
field is calculated to be 2.6 G, based on the coil geom- 
etry. We neglect the resulting differential Zeeman shifts 
of 600 kHz [19 because they are small compared to the 



cavity linewidth of 5.3 MHz. 

The experiments were performed on the transitions 
from the \F = 2) ground state to the \F' = 1,2,3) ex- 
cited states of the D2 line [19 . A model Hamiltonian 
describing the system is 

N 3 

R=huJc{a^a^b^b)^hY^ ^ \F') 



i=l F' = l 



N 3 



E {9^a^D^i2,F')+g,b^Dli2,F')+h.c.). 

(1) 



i=l F' = l 



Here, a and b are the annihilation operators for the cavity 
modes corresponding to tt and ± transitions of the atom, 
ujc is the resonance frequency of the cavity, \F') {F'\. is 
the operator projecting the i-th atom onto the excited 
state \F')^ uoy' is the frequency of the transition from the 
\F = 2) ground state to the \F') excited state, gi is the 
complex atom-cavity coupling constant of the z-th atom 
and h.c. denotes the Hermitian conjugate. D^{2^F') and 
D^{2^F') are the atomic dipole transition operators for 
the z-th atom from the \F') excited state to the \F = 2) 
ground state, which take into account the different cou- 
pling strengths for the tt and ± transitions. 

We follow the convention in [4 to define the atomic 
dipole transition operators for an atom interacting with 
different polarizations of the light field as 

D\F,F') = 

|F, tuf) {F, mF\fiq\F\ ruF + q) (F^ + g| , (2) 

where q = { — 1,0,1} and /ig is the dipole operator 
for TT, (j+}-polarization, normalized such that for 

the cycling transition from the \F = 2,mi? = 2) to the 
\F^ = 3, mi?/ = 3) state (/i) = 1. For ^-polarized light 
we identify = {D^^ + A"^)/V^- 

In the experiment, we probe the system with tt- 
polarized light after having prepared the atoms in a sta- 
tistical mixture of mi? states. For weak driving and large 
atom numbers the system behaves like a collective spin 
with one ground and three excited states coupled to the 
TT-mode of the cavity, as indicated in Fig.[l](b). We model 
the system by the effective Hamiltonian 



H =nujca^a -^Ti^uj^' \F') {F'\ 



F' = l 



n (yNgp'a^ \F = 2) + h.c.) , 



(3) 



which we derive in detail in Appendix [A| The coupling 
of the collective spin is enhanced by a factor as com- 
pared to an effective single atom coupling constant gF' 
that is an average of both the spatial dependence and 
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the rriF state dependence of the couphng of the indi- 
vidual atoms. Because we are only considering the case 
of weak driving, we restrict our discussion of the spec- 
trum to the excitation manifold with a single excitation 
present in the system [4 . In this manifold, there are three 
avoided crossings between the energies of the bare cavity 
and the bare atom. If the crossings are well separated, 
we expect three splittings in the transmission spectrum 
of size 2y/NgF'^ separated by the hyperfine splitting of 
the excited states of ^'^Rb. Solving the eigensystem for 
the Hamiltonian yields a more detailed description of the 
spectrum. We expect to see transmission when the probe 
laser frequency is resonant with the eigenenergies of the 
system, while the amount of transmission is proportional 
to the expectation value {a^ a) of the corresponding eigen- 
state [J. Numerical calculations lead us to predict the 
detailed form of spectrum as summarized in Fig. [2] 
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III. EXPERIMENT 

We verify our predictions using the experimental setup 
sketched in Fig.[l](a). Its central element is a cavity made 
of two coned down spherical mirrors with radii of curva- 
ture of 25 mm separated by 500 /im. The transmission 
per mirror is ~ 20 ppm and the losses are ~ 35 ppm, 
as estimated from measurements of the transmission, re- 
flection and linewidth of the cavity. In terms of cavity 
QED parameters, our setup is described by (^, hz^ 7) = 
27r X (9.2,2.6,3.0) MHz, where g is the maximum sin- 
gle atom coupling constant for the \F = 2,mF = 2) to 
\F^ = 3, mi? = 3) transition, is the cavity field decay 
rate, and 7 is the atomic dipole decay rate p!9 | [2Q j . 

We start the experiment by loading atoms into an op- 
tical dipole trap in the |F = 1) ground state, similar to 
our previous experiments [21 . The dipole trap is formed 
by a beam with a wavelength of 1064 nm and a power of 
12 W focused to a waist of 25 /im. From there the atoms 
are evaporatively cooled into a transfer lattice formed by 
two counter-propagating beams of the same wavelength. 
Each lattice beam has a maximum power of 1 W and is 
focused to a waist of 50 /im inside the cavity. By chang- 
ing their relative frequency, the atoms are transported 
into the cavity [22|, which is located 9.2 mm below the 
position of the dipole trap. Inside the cavity, the atoms 
are trapped in a two-dimensional lattice formed by the 
transfer lattice and an intra-cavity beam at 808 nm with 
a circulating power of 3 mW and a waist of 25 jam. We 
also use this beam to stabilize the length of the cavity. 

Within the cavity, the atoms are repumped into the 
\F = 2) ground state and the repumping beam is left on 
during the rest of the experiment. The atom-cavity sys- 
tem is probed through the cavity with light linearly po- 
larized along the axis of quantization defined by the mag- 
netic field. The intensity of the probe is adjusted to give 
a small mean photon number nempty in the empty cavity. 
It is chosen to be as close as possible to the weak prob- 
ing condition nempty <^ 1 while still resulting in enough 



FIG. 3. Cavity transmission. For a given experiment we 
fix the cavity detuning Ac and the atom number N. The 
probe laser detuning is swept from —700 MHz to +700 MHz 
in 210 ms, which is small compared to the measured atom 
number lifetime of 3 s. At the end of the sweep the atom 
number is measured by absorption imaging. The cavity out- 
put is directed onto a single photon counting module and the 
counts are binned in 100 /is intervals. All detunings are with 
respect to the |F = 2) to = 3) transition. 



signal for the experiment. The transmission spectrum is 
recorded by directing the output of the cavity onto a sin- 
gle photon counting module. After probing the system, 
we measure the atom number via absorption imaging. 

For each experiment the cavity detuning Ac = uq—^i 
is fixed and we control the number of atoms by adjusting 
the power in the initial dipole trap. During the exper- 
iment the probe laser frequency is swept over the fre- 
quency range of interest, resulting in a transmission sig- 
nal as shown in Fig. [3j By repeating the experiment 
for different cavity detunings and different average atom 
numbers, we map out the transmission spectra, as shown 
in Fig. [4) 



IV. CONCLUSION 

To compare the measured spectra with our model pre- 
dictions, we assume equal population of both the mp 
states and the lattice sites of the intra-cavity trap. For 
simplicity, we neglect the finite size of the cloud along 
the direction transverse to the optical axis of the cavity 
[23^. Because the intra-cavity trap has a different lattice 
spacing than the mode of the cavity, the spatial aver- 
aging of Qi gives g ~ m8ix{gi) / \^ = 6.5 MHz, as for a 
uniform distribution of atoms along the cavity axis. The 
agreement between the experiment and the model, as il- 
lustrated in Fig. [2] and Fig. [4j is obtained without any 
free parameters. 
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FIG. 4. Experimental transmission spectra. By varying 
the cavity detuning we map out the full transmission spec- 
tra of the system for different average atom numbers N, 
with their standard deviation indicated in parentheses. The 
recorded transmission is indicated in color in terms of intra- 
cavity photon number n, normalized to the average photon 
number inside the empty cavity on resonance riempty, where 
^empty = 0.25 for the top left and bottom right plots and 
^empty = 0.35 for the top right and bottom left plots. Points 
with counts below 2 are colored as 0. The gray solid lines 
indicate the positions of the theoretical eigenenergies for the 
measured average atom number an equal distributions of 
TTiF states and a ^ of 6.5 MHz. All detunings are with respect 
to the |F = 2) to |F' = 3) transition. 



So far, we have neglected the effect of non-neglible driv- 
ing strength. Our model assumes weak driving in order 
to describe the system as a collective spin. It will fail 
for higher driving strength since it predicts a strong non- 
linearity due to the saturation of the single spin. How- 
ever, as the atom number increases, the system saturates 
at higher excitation numbers and loses the strong non- 
linearity predicted by our simple model. Another effect 
we have neglected are the tuf state changing processes 
either due to optical pumping or due to the weak re- 
laxation via the undriven mode of the cavity. This is 
justified because the experiment was performed with a 
large number of atoms in a statistical mixture of all rriF 
states. Thus, residual tuf state changing processes are 
not expected to alter the distribution significantly. In a 
series of related experiments, we have started to inves- 
tigate the dynamics of polarized gases, where the atoms 
are optically pumped into a particular ttif state prior to 
their interaction with the cavity. In this case, we see ev- 
idence of mi? changing processes, which in turn change 
the effective coupling qf' of the system. These processes 
lead to complex dynamics, which will be an interesting 
topic for future investigations. 

In summary, we have experimentally demonstrated 
that the linear response of a gas of alkali-metal atoms 



coupled to a high-finesse optical cavity is well described 
by a collective spin with multiple levels coupled to a sin- 
gle mode of the cavity. Using this system provides a 
flexible testbed for collective cavity QED with multiple 
atomic levels and the theoretical frame work is appli- 
cable to other physical systems used in quantum infor- 
mation processing. These include hybrid systems like 
alkali-metal atoms above a coplanar waveguide, collec- 
tively coupled to linearly polarized microwave photons 
[24] . and nitrogen- vacancy centers coupled to a super- 
conducting resonator, where a similar effect in the trans- 
mission spectrum has recently been observed [25] . 
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Appendix A: Derivation of the reduced Hamiltonian 

In order to derive the effective Hamiltonian of Eq. [3j 
we take a closer look at the interaction part Hi of the full 
Hamiltonian of Eq. [l] that is, 

N 3 

Hi = E E (5i«^A°(2,F') +5^6tA^(2,F') +h.c.) . 

i=l F' = l 

(Al) 

We will first consider the driven 7r-mode and then the un- 
driven ±-mode. Driving the 7r-mode of the cavity pop- 
ulates the state where all the atoms are in the ground 
state and there is one 7r-photon in the cavity. 



N 



k,gAr) = k) JJ|F,mF)i 



(A2) 



Subsequent interaction with the cavity leads to transi- 
tions to the Dicke state, where one of the atoms is excited 
and the cavity is empty, 

|0,e^) = |0) |^^El^''"^^')«ni^'"^^)ij ■ (A3) 



The rate of which is collectively enhanced, i.e., 
(0,eAr|Hi|7r,gAr) = \fNgF'' 



(A4) 



To calculate the effective single atom coupling rate qf' , 
we note that the coupling of the atom to the 7r-mode of 
the cavity depends on its position relative to the cavity 
field. Along the cavity axis the corresponding Qi varies 
from maximum coupling for an atom at an anti-node of 
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the cavity field, to zero coupling for an atom at a node. In 
the transverse direction the coupling varies as the Gaus- 
sian field distribution of the TEMqo mode of the Fabry- 
Perot cavity. 

In addition to position, the coupling depends on the 
rriF state of the atom via the Clebsch-Gordon coefficient 
(F, mi?|/io|^^ ^f)^ that enters the dipole transition op- 
erator 

dUf,f') = 

^ |F, {F, mpltiolF', mF)^ {F\ mF\i ■ 

rriF 

For a product state of N atoms, each in a particular tuf 
state, the sum in the dipole transition operator D^{2^ F')i 
has only one non-zero term for every atom. The oper- 
ator then becomes a lowering operator, taking the i-ih 
atom from the \F') to the \F = 2) state, with a prefactor 
given by the Clebsch-Gordon coefficient and the spatially 
varying coupling constant gi^ i.e., 

giD°{2,F')i = gi (F, |mo|F', \F = 2) . 

(A6) 

For weak driving the ensemble behaves as a collective 
spin and the effective single atom coupling constant re- 
sults from both an average over the spatial distribution 
of the atoms and their tuf states, that is 



\ 



1 ^ 

-Y,\{F.mF\^o\F\mF),g^\' (A7) 



N 



'gJ^p{mF)mmF\fio\F^,mFW. (A8) 



population in the different ttif states. In the last step, 
we have assumed that enough atoms are in all of the ttif 
states such that the atoms in each state independently 
average to g. 

The undriven mode of the cavity is only populated 
by a transition of the excited state |0,eAr) into the state 
where the cavity holds a ±-photon and the atoms are in a 
superposition of one of them having changed tuf states. 



I.e., 




F,mF^q),Y[\F,mF)^ . 



(A9) 

The rate at which this process occurs is not collectively 
enhanced [l6l[T8], i.e.. 



(^,giv|Hi|0,eAr) = ^f', 



(AlO) 



even if all tuf states are macroscopically occupied jTT] . 
Because ^ 1 in the experiment, we restrict our treat- 
ment to the interaction with the driven 7r-mode of the 
cavity and arrive at the effective Hamiltonian for weak 
driving 



H =hujca'^a + ^ cjp' \F') {F'\ 



F' = l 



h iyNgp'o) \F = 2) (F'l + h.c.) . 



(All) 



F' = l 



Here, g = \/Xlil^*P/A^ the coupling constant resulting 
from spatial variations in gi and p{mF) is the relative 



The treatment above is similiar to work in the area of 
quantum dots studying the vacuum-Rabi splitting in the 
presence of inhomogenous broadening [26] . 
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